Linear-response theory of the longitudinal spin Seebeck effect 
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We theoretically investigate the longitudinal spin Seebeck effect, in which the spin current is 
injected from a ferromagnet into an attached nonmagnetic metal in a direction parallel to the 
temperature gradient. Using the fact that the phonon heat current flows intensely into the attached 
nonmagnetic metal in this particular configuration, we show that the sign of the spin injection 
signal in the longitudinal spin Seebeck effect can be opposite to that in the conventional transverse 
spin Seebeck effect when the electron-phonon interaction in the nonmagnetic metal is sufficiently 
large. Our linear-response approach can explain the sign reversal of the spin injection signal recently 
observed in the longitudinal spin Seebeck effect. 
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Because of the desire to deal with heating problems 
in modern spintronic devices, there has been an increas- 
ing interest in investigating thermal effects in spintron- 
ics. A new subfield "spin caloritronics" [l[ aims to un- 
derstand the basic physics behind the interplay of spin 
and heat. One of the central issues in spin caloritronics 
is the newly discovered thermo-spin phenomenon termed 
spin Seebeck effect Q, which enables the thermal injec- 
tion of spin currents from a ferromagnet into attached 
nonmagnetic metals over a macroscopic scale of several 
millimeters. The spin Seebeck effect is now established 
as a universal aspect of ferromagnets because this phe- 
nomenon is observed in various materials ranging from 
the metallic ferromagnets Ni 81 Feig @ and Co 2 MnSi Q, 
to the semiconducting ferromagnet (Ga,Mn)As [J], to the 
insulating magnets LaY2Fe50i2 [5(. 

It is important to note that the above experiments 0- 
Q were performed in a configuration of the transverse 
spin Seebeck effect, in which the direction of the ther- 
mal spin injection into the attached nonmagnetic metal 
is perpendicular to the temperature gradient [Fig. [T] (a)]. 
Recently, another type of spin Seebeck effect called the 
longitudinal spin Seebeck effect 0,0) is reported, in which 
the direction of the thermal spin injection into the non- 
magnetic metal is parallel to the temperature gradient 
[Fig. 0(b)]. Whereas the longitudinal spin Seebeck effect 
is well defined only for the use of an insulating ferromag- 
net due to theparasitic contribution from the anomalous 
Nernst effect [1, @, it has several attractive features: (i) 
it is substrate free, (ii) the configuration is much simpler 
than that of the transverse spin Seebeck effect, and (iii) 
it can be of wide application because it allows the use of 
bulk samples. 

Another pronounced feature of the longitudinal spin 
Seebeck effect is that the sign of the spin injection sig- 
nal is opposite to that in the transverse spin Seebeck cf- 
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FIG. 1: (Color online) Schematic view of the experimental 
setup for (a) the transverse spin Seebeck effect and (b) the 
longitudinal spin Seebeck effect. 



feet [f| 0] ■ Physically, the longitidinal spin Seebeck effect 
is distinguished from the transverse spin Seebeck effect by 
the fact that the attached nonmagnetic metal is in con- 
tact with the heat bath in the longitudinal setup, while 
the attached nonmagnetic metal is out of contact with 
the heat bath in the transverse setup. This brings about 
a clear difference that the heat current intensely flows 
into the attached nonmagnetic metal in the case of the 
longitudinal spin Seebeck effect, whereas it does not in 
the case of the transverse spin Seebeck effect. It is obvi- 
ous that theory of magnon-driven spin Seebeck effect 
fails to explain the situation in question. 

In this paper, by em ploy ing linear-response theory of 
the spin Seebeck effect [Tljj and using the importance of 
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the phonon-drag process in the spin Seebeck effect [12J, 
we show that the sign of the spin injection signal in the 
longitudinal spin Seebeck effect can be opposite to that in 
the conventional transverse spin Seebeck effect when the 
electron-phonon interaction in the attached nonmagnetic 
metal is sufficiently large. The key in our discussion is the 
aforementioned difference in the position of the attached 
nonmagnetic metal between the longitudinal setup and 
the transverse setup. 



II. PHENOMENOLOGY OF THE 
LONGITUDINAL SPIN SEEBECK EFFECT 

Let us begin with the phenomenology of the longitu- 
dinal spin Seebeck effect. In Fig. [TJ a hybrid structure 
of a ferromagnet (F) and a nonmagnetic metal (A) is 
placed under a temperature gradient. The central quan- 
tity that characterizes the spin Seebeck effect is the spin 
current I s injected into N. As explained in detail in 
Ref. the spin Seebeck effect is a thermal spin injec- 
tion by localized spins, and the injected spin current has 
two contributions, 
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where /pump ^he so-called pumping component) repre- 
sents the spin current pumped into A by the thermal fluc- 
tuations of localized spins in F , while i s back (the so-called 
backffow component) represents the spin current coming 
back into F by the thermal fluctuations of the spin ac- 
cumulations in A. We now focus on the spin current in- 
jected into A which is located close to the cold reservoir. 
In the case of the conventional transverse spin Seebeck 
effect, the pumping component JP um P is lesser than the 
backffow component 7 s pump [Fig-COa)]. In contrast, /P um P 
is greater than the backffow component /pump m the case 
of the longitudinal spin Seebeck effect [Fig. O^b)]. 

This difference can be explained phenomenologically 
on the basis of the following conditions: (i) most of the 
heat current in the F/N hybrid system at room tem- 
perature is carried by phonons (see Ref. [14| in the case 
of yttrium iron garnet), and (ii) the interaction between 
the phonons and the spin accumulation in A is much 
stronger than the magnon-phonon interaction in F. 

First, recall that the pumping and backflow compo- 
nents can be expressed as follows [l3| : 
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where T F and are the effective temperature of the 
magnon in F and the spin accumulation in N. Here, 
Gs = JsdXNT s {/H with J sd , xn, and r sf being the s-d in- 
teraction at the interface, the paramagnetic susceptibil- 
ity in A, and the spin diffusion length in A, respectively. 
The negative sign before G s arises from the fact that the 
magnon carries spin — 1. In the longitudinal spin See- 
beck experiment, the nonmagnetic metal A is in direct 



contact with the heat bath, and thereby is exposed to 
the flow of the phonon heat current due to condition (i). 
Then, because of condition (ii), spin accumulation in A 
is heated up faster than the magnons in the ferromagnet 
F, and the resultant effective temperature of the spin 
accumulation in A increases above that of the magnons 
in F. In the conventional transverse spin Seebeck setup, 
by contrast, the nonmagnetic metal A is out of contact 
with the heat bath and the phonon heat current does not 
flow through the nonmagnetic metal A, while the ferro- 
magnet F is in contact with the heat bath, resulting in 
an increase in the effective magnon temperature in F. 
Therefore, in this case, the effective temperature of the 
spin accumulation in A is lower than that of the magnons 
in F. This difference can explain the sign reversal of the 
spin Seebeck effect signal between the longitudinal setup 
and the conventional transverse setup. 



III. LINEAR-RESPONSE FORMULATION 

In this section we review the linear-response formalism 
of the spin Seebeck effect developed in Ref. [ll|- In the 
next section, this formalism is employed to evaluate the 
longitudinal spin Seebeck effect. We use a model shown 
in Fig. [21 in which the localized spins in F are interact- 
ing with the spin accumulation in A through the s-d ex- 
change interaction J sc j at the interface. In our approach, 
the spin accumulation is modeled as a nonequilibrium 
itinerant spin density s. 

As in Ref. [ill ], the spin current I s injected into 
the nonmagnetic metal A is calculated as I s (t) = 
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number of lattice sites in F (A), So is the size of the 
localized spins in F, and l7 S( j +9 is the Fourier trans- 
form of the s-d interaction at the F/N interface. Here, 
C£ (t,t') = — i(a+(i')Sk (t)) measures the correlation be- 
tween the magnon operator a+ in F and the itinerant 
spin-density operator s^ = (sj£ — is|)/2 m N. Note that 
the time dependence of I s (t) vanishes in the steady state 
and it is hereafter discarded. Introducing the frequency 



representation C< q (t - f ) = %C< h (u) 



and adopting the representation [l5| C 

well as using the relation C < — >ir ' k 
obtain 
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for the spin current I s in the steady state [l6| . 

Up to the lowest order in the s-d interaction J s( j, the 
interface correlation function C appearing in Eq. (|4]) is 
generally expressed as 
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FIG. 2: (Color online) Diagrammatic representation of the 
spin Seebeck effect in the longitudinal configuration. The 
double solid line (the dashed line) represents the propagator of 
the itinerant spin density (the phonon) . The solid circle (solid 
triangle) denotes the interaction vertex TK, q (T^.fc) between 
the magnon and the phonon (the itinerant spin density and 
the phonon). Here, T\ < T 2 < T 3 . 



where X q (uj) — X q (uj) + 5X q (uj) is the renormalized 
magnon propagator with the bare component X q (u), and 
X fc (w) = Xk(w) + Sxk(w) is the renormalized spin-density 
propagator with the bare component Xfc( w )- The bare 
magnon propagator satisfies the equilibrium condition: 



X£(u>) = [X?(lj)}*, = 2i Ii<( W ) coth( 



(6) 



where the retarded component is given by X q (u) = 
(uj — uiq + icoj) -1 with uj q = 7-ffo + w q being the magnon 
frequency for uniform mode jHq and exchange mode 
uj q = D cx q 2 /h. Likewise, the bare spin-density propa- 
gator satisfies the local equilibrium condition: 

X* M = [X?M]*, X* M = 2iIm X ^(c) coth(^), 

(7) 

where the retarded component of Xfc( w ) i s given by 
Xfc (w) = Xjv/(1 + ^sffc 2 — iwrrf) with A s f being the spin 
diffusion length. 



IV. CALCULATION OF THE LONGITUDINAL 
SPIN SEEBECK EFFECT 

In this section we present a linear-response calculation 
of the longitudinal spin Seebeck effect and justify the 
phenomenological picture presented in Sec. [TTJ The spin 
current I s injected into N due to the longitudinal spin 
Seebeck effect is composed of three terms, 



7.(o) +/,(&)+ 7. (c), 



(8) 



where I s (a), 7 S (6), and 7 s (c) correspond to diagram (a), 
(b), and (c) in Fig.[5J Below we show that each term has 
the following sign: 



la (a) < 0, I s (b) < 0, 7 s (c) > 0. 



(9) 



First, let us consider diagram (a). This can be calcu- 
lated by setting ^(w) — ► X q {uj) and X fc ( w ) — > Xfc( w ) 
in Eq. (J5)), which was already done in Ref. [ll| with the 
injected spin current given by (see Eq. (12) therein) 
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where we have introduced the shorthand notation f = 

J u 

2 ' ana - "^ int * s * ne num b er °f localized spins at the 
N/F interface. Note that I s (a) has a negative value due 
to T 2 > Ti. 

Next, let us consider diagram (b). In this process, 
the localized spins in F is excited by the nonequilibrium 
phonon driven by the temperature gradient in F, hence 
this corresponds to the phonon-drag process [l2j. Eval- 
uation of the diagram was already given in Ref. [Hj], al- 
though the calculation is lengthy and tedious (see the 
supplemental material therein). In short, this term can 
be calculated by setting in Eq. © X q {uj) -> SX'^Iuj) = 

(q' ) with its Keldysh component given by 
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8X^' k (uj) = -2V-^ / 5D 
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x |^H| 2 [coth(^)]-coth(^) 

(11) 

where we have introduced shorthand notations cj_ = oj — 
u,q- = q-K, and /„ = fg, and 



K,q — 9n 



2M ion ^ h 



(12) 



is the magnon-phonon interaction vertex. Here, vk, 
v p h and M- lon are the phonon frequency, phonon ve- 
locity and the ion mass, respectively, and the strength 
of the magnon-phonon coupling is given by g m -p = 
|asV Jex\(wq/ J e x) with the exchange interaction J cx . In 
Eq. HID, 

5D n K eq ' K {v) = 2i 



\n K+K '\ 2 



K' 



x [eoth(^)-coth(^)] (13) 

is the Keldysh component of the nonequilibrium phonon 
propagator. Here Q^~ K is the Fourier transform of 

ft ph (r) = ^0Er o eF/ J F-intcrfacc4 <5 ( r ~ r o) With ^0 = 
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\/2ifph/Mion, K P h is the elastic constant, and a| is the 
cell volume of the ferromagnet. Putting these expressions 
into Eq. ([5]) and after some algebra, we finally obtain 
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F k,q,K,K' Jv 

xImZ?£>)[coth(^) -coth(^)^4) 



where -Dg-(^) = (^-^ + i/r ph ) 1 - (V + ^jt + i/r ph ) Ms 
the retarded component of the phonon propagator with 
the phonon lifetime r ph , X = ^/2(J^ d S )CllN int N[ n jN F 
with A^ nt being the number of lattice sites at the F/F 
interface, and Ak, q {v) is defined by 

A k , q {v) = f lm X %[w)lmX*_{u-) 

x |X«H| 2 [coth(^)-coth(^)]. (15) 

Note that only the even component of Ak jq (v) as a func- 
tion of v gives a non- vanishing contribution to Eq. (|14p. 
Because the even component of Ak, q {v) is negative def- 
inite as well as lm.D^(u) [coth( 2 ^ Ta ) - coth( ^^r )] in 
Eq. (TH)) . I s (b) has a negative value. 

Finally, let us consider diagram (c). Repeating essen- 
tially the same procedure in evaluating diagram (b), we 
obtain 
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E ( T ^) 2 / B ktq {v)\D*(y) 



xIm^«,H[coth(^) - coth(^)]Gl6) 

where V = y/2{J^ d S Q )£llN? nt /N F , D%{v) denotes the 
phonon propagator in N . In the above equation, the cou- 
pling between the itinerant spin density and the phonon 
in N is given by 
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2M; C 



2 < 



(17) 



where g s - p s» |aVtho P |C^ 2 ^V'(0) with a and ihop being the 
lattice spacing and the hopping integral of the nonmag- 
netic metal N. In Eq. (flTJl) . Bk, q {v) is defined by 

x[coth(^)-coth(^)]. (18) 

Note that as in Eq. (fT4|) . only the even-in-z/ component 
of Bk, q {v) gives non- vanishing contribution to Eq. (1161) . 
Then, because the even-in-^ component of Bk. q iy) is neg- 
ative definite, I s {c) has a positive value. 

V. DISCUSSION 

In the previous section, we proved that I s {a) and I s (b) 
have the same sign, whereas I s (c) have the opposite sign 



[Eq. ([§)]. Then, if I s (c) is dominant in Eq. ((5J, it means 
that the sign of the longitudinal spin Seebeck effect can 
be opposite to that of the transverse spin Seebeck effect, 
since the sign of I s (a) and I s (b) is the same as that of 
the transverse spin Seebeck effect (see Refs. [ll| and |12|). 
Because I s (a) and I s (b) are considered to have the same 
magnitude at room temperature (see Fig. 3 in Ref. 
we here compare the magnitude of I s (b) and I s {c). 

The key quantities determining the magnitude of 
I s (b) and I s (c) are the interaction vertex Tjc,q between 
magnons and phonons [solid circle in Fig. [2] and Eq. (|12l) ] 
and the interaction vertex T k ,k between spin accumula- 
tion and phonons [solid triangle in Fig. [2] and Eq. (|17)) ], 
The magnitude of these couplings is roughly given by 
9m- P » |asVJ ex |(wq/Jex) and g s - p fa |aVi h op|C/ 2 -/V'(0), 
where a, ihop, U, N'(0) are the lattice spacing, the hop- 
ping integral, the strength of the Coulomb repulsion, 
and the strength of the particle-hole asymmetry in the 
nonmagnetic metal N. For materials with a relatively 
large Coulomb repulsion U [H[ and moderate strength 
of particle-hole asymmetry N'(0) [l9| such as Pt, we ex- 



pect a situation g s 



> 



im—pt 



which then explains the 



sign reversal of the spin injection signal in the longitu- 
dinal spin Seebeck effect due to Eqs. dHJ and ((9|). From 
these considerations, we conclude that this happens in 
the longitudinal spin Seebeck effect reported in Refs. @ 
and 01. 



VI. CONCLUSION 

In this paper we have developed linear-response theory 
of the longitudinal spin Seebeck effect. We have shown 
that the sign of the spin injection signal in the longitu- 
dinal spin Seebeck effect can be opposite to that in the 
conventional transverse spin Seebeck effect when the in- 
teraction between the spin accumulation and the phonon 
in the attached nonmagnetic metal is sufficiently stronger 
than the interaction between the magnon and the phonon 
in the ferromagnet. The linear- response approach pre- 
sented in this paper can explain the sign reversal of the 
spin injection signal recently observed in the longitudinal 
spin Seebeck effect [f| 0] • 
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Appendix A: Calculation of the vertex T 

In this Appendix, we evaluate the interaction vertex 
T K,k (solid triangle in Fig. [5]) between the itinerant spin 
density and the phonon in the nonmagnetic metal N . We 
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FIG. 3: Diagram corresponding to the interaction vertex 
Tic.te between the itinerant spin density and the phonon. 
The solid line, the wavy line, the dashed line and the triple 
dashed line represent the electron propagator, the Coulomb 
repulsion, the phonon propagator, and the diffuson propaga- 
tor, respectively. "R" or "A" means the retarded or advanced 
component of the electron Green's function. 



assume that the nonmagnetic metal N has a moderately 
large Stoner enhancement factor, and for conduction elec- 
trons in N we use a model described by Ref. and 
assume an elastic impurity scattering as well. The inter- 
action vertex before integrating out the fermionic degrees 



of freedom is shown in Fig. [3l The building block of this 
diagram is given by a triangle 



T = 



tanh( 



2k B T 



I lmG*_ h (e- U )G*{e)G$(e),(AL) 

J n 



where Gp^ A (e) = (e — e p ± i/r) -1 is the electron Green's 
function with the electron's lifetime r. This diagram can 
be evaluated to be T ~ N'(0)ujt as was done in Ref. [20l |. 
After the inclusion of a diffuson vertex correction (triple 
dashed ladder in Fig. [3]) which is important in a realistic 
diffusive situation, we obtain 



T w JV'(O) 



{Dk 2 r) 2 + {wrf 



(A2) 



where D is the diffusion constant. The dominant contri- 
bution comes from the dynamical region u 3> Dk 2 and in 
this case we approximately have T ~ iV'(O). By attach- 
ing two Coulomb repulsion U and one electron-phonon 

interaction w |aVfhopL / 2m" K v' 2 cormn g from each 

' "* ph 



vertex, we finally obtain Eq. (fT71) . 
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